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Phase separation in a complex fluid with lamellar order has been studied in the 
case of cold thermal fronts propagating diffusively from external walls. The veloc- 
ity hydrodynamic modes are taken into account by coupling the convection-diffusion 
equation for the order parameter to a generalised Navier-Stokes equation. The dy- 
namical equations are simulated by implementing a hybrid method based on a 
lattice Boltzmann algorithm coupled to finite difference schemes. Simulations show 
that the ordering process occurs with morphologies depending on the speed of the 
thermal fronts or, equivalent ly, on the value of the thermal conductivity ^. At large 
value of ^, as in instantaneous quenching, the system is frozen in entangled config- 
urations at high viscosity while consists of grains with well ordered lamellae at low 
viscosity. By decreasing the value of ^, a regime with very ordered lamellae parallel 
to the thermal fronts is found. At very low values of ^ the preferred orientation is 
perpendicular to the walls in d = 2, while perpendicular order is lost moving far 
from the walls in d = 3. 
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1. Introduction 

Ordering in complex fluids can be a slow process. Thermodynamic equilibrium is 
often characterised by the presence of organised structures at mesoscopic scales 
(Larson 1999). However, the dynamical formation of topological defects can make 
difficult to reach the true equilibrium state (Boyer & Viiials 2002, Benzi et al. 2009). 
Since there are physically relevant properties related to the degree of order in the 
system, it can be crucial to control the ordering process by proper driving. 

In this paper we will consider fluids with lamellar order, such as symmetric 
block copolymer melts (Fredrickson & Bates 1996). For recent applications of these 
systems in micro- and nano-technologies one can for example see the work by Man 
et al. (2010) where the relevance of having long-range order is discussed. A tran- 
sient macroscopic sample of a copolymer melt usually consists of different oriented 
lamellar domains (grains) with defects like dislocations and disclinations. The in- 
terplay between the mesoscopic structures and the local velocity field is essential 
for the evolution of defects (Gonnella et al. 1997). Hydrodynamic modes favour the 
disentanglement of intertwined patterns but simulations have shown how extended 
defects between different grains remain at late times slowing down the dynamics of 
ordering (Xu et al. 2004, 2005). 
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Mechanical driving like shear can favour ordering since interfaces generally pre- 
fer to be aligned with the flow (Larson 1999). In lamellar systems under shear, 
however, different stable orientations arc possible (Frcdrickson 1994) and in the 
ordering process (Gates & Milner 1989, Corberi et al. 2002) the shear by itself can 
cause the creation of defects due to the dilatation or compression of the layers (Ku- 
maran et al. 2001, Lu et al. 2008). More complex phenomena like shear banding 
have been also observed (Xu et al. 2006). 

A different controlling procedure could be thermal driving. Recently, in simple 
binary mixtures, it has been shown than in thermally controlled quenches with cold 
fronts propagating into the system at finite speed, interfaces have a preferred orien- 
tation with respect to the fronts (Furukawa 1992; Foard & Wagner 2009; Gonnella 
et al. 2010). As a result, ordered patterns arc obtained. Applications to polymer 
systems have been considered for example by Voit et al. (2005). In this work we 
want to investigate how a thermally controlled quench could be used to produce 
well ordered lamellar systems. 

We consider a system described in the continuum by the Navier-Stokes and 
convection-diffusion equations with the thermodynamics of the lamellar phase en- 
coded in a free-energy functional. The temperature evolves by following a standard 
diffusion equation with the system between two walls at a temperature below the 
critical value. In order to solve the dynamical equations we apply a hybrid method 
with a lattice Boltzmann algorithm for solving the Navier-Stokes equations and a 
finite difference scheme for the convection-diffusion equation. The model and the 
method will be described in the next section while the numerical results will be 
presented Section III. Few conclusions will complete the paper. 



2. The model 

We consider a binary mixture having total density n and density difference (p, 
the order parameter, between the two components. The equilibrium properties are 
described by the free energy (Brazovskii 1975) 



dr. (2.1) 



The first two terms depending on n and on the temperature T (the Boltzmann 
constant is assumed to be one) give the ideal gas pressure which does not affect the 
phase behaviour. We consider 6, c > to ensure stability. The parameter a can be 
thought of as a reduced temperature 

T — T 

a = i— ^ (2.2) 

c 

and is the only one depending on temperature in the terms depending on ip. For 
temperatures above the critical value the fluid is disordered. For T < and 
K > there is the coexistence of two homogeneous phases with equilibrium val- 
ues if = —a/b. Negative values of k favour the formation of interfaces, while 
a reduction of n can induce a transition into a lamellar phase. By using a single 
mode approximation with proflles like A sin Uqx in the direction normal to lamel- 
lae, one flnds a transition, when \a\ = 6, at a ~ — I.IIk^/c with fco = \J —Kjlc and 
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= 4(1 + K^/4c6)/3 (Xu et al. 2005). The parameters of this model have been 
mapped on those of copolymer systems by Binder (1995). 

The evolution of the system is described by a set of coupled partial differential 
equations: 

dtn + da(nua) — 0, 



n{dtUa + updpUa) 



dt(p + da{(pUa) = rV^^, 



(C - 2ri/d)dai3djU^ 



(2.3) 
(2.4) 
(2.5) 



where u is the local velocity of the fluid, 77 the shear viscosity, C the bulk viscosity, 
d the space dimension, and F the mobility coefficient. We model the quench process 
by taking the reduced temperature a (j2.2[) to be the solution of a diffusion equation 
dta — ^V'^a where ^ is the thermal conductivity. The pressure tensor Pais and the 
chemical potential fj, can be computed from the free energy functional (j2.ip and 
have the following forms (Gonnella et al. 1997): 



Pa0 ^ [nT + -^^ + -^^ 



^(VV) + i(v^)2 



■^(v^)V + ^(vV)^ 



+djipdj{V'^ip) ^Sap + Kdaifdpif - c daipdfj{V^(p) + dfj(pda{V^ip) (2.6) 



and 



(2.7) 



Equations (|2.3p - (|2.5p are solved numerically, by using a mixed approach. A lat- 
tice Boltzmann scheme is used for the continuity and Navier-Stokes equations (I2.3P 
and (|2.4p . and a finite-difference approach for equation (|2.5I) and for the temper- 
ature diffusion equation. The full study is limited to bi-dimensional systems while 
the equations for ip not coupled to hydrodynamics will be solved also in cZ = 3. 
In the case of a simple fluid (Benzi et al. 1992, Chen & Doolen 1998, Succi 2001) 
the lattice Boltzmann method (LBM) is defined in terms of a set of distribution 
functions, fi{r,t), located in each lattice site r at each time t, and of a set of nine 
velocity vectors e^, defined on a square lattice, having moduli |ei| = 0,c, V^c with 
c = Ax/ At being Aa; and At the lattice and the time step, respectively. The distri- 
bution functions evolve according to a single relaxation time Boltzmann equation 
(Bhatnagar et al. 1954) 



Mr + e,At,t + At)~ f,{r,t) 



At 



[h{r,t)~ft\r,t)]+AtF.,{v,t), (2.8) 



where r is the relaxation time, f^'' are the equilibrium distribution functions and 
Fi are the forcing terms to be properly fixed. 

The physical quantities, that is the total density n and the fluid momentum nu, 
are determined by the relations 



nu = ^ /,e, + ipAt, 



(2.9) 



with F the force density acting on the fluid. The equilibrium distribution functions 
f^'^ are given by the standard second order expansion in the fluid velocity u of the 
Maxwell-Boltzmann distribution function (Qian et al. 1992). 
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The forcing term in equation (I2.8P is expressed as a second order expansion in 
the lattice vector velocities (Guo et al. 2002) 



F,,= 



Gj U Gi • U 

^ ^ e,; 



•F, (2.10) 



where Cs = c/a/3 is the speed of sound and uji — 4/9, 1/9, 1/36 in correspondence 
of the lattice directions with |ei| — 0, c, \/2c, respectively. The force F has to have 
the following expression 

Fp = d^{nc%p - P^p) (2.11) 
in order to recover equation (j2.4l) where it results 

C - . - nciAt - 1) . (2.12) 

The first term on the r.h.s. of equation (|2.1ip allows to cancel out the athermal ideal 
gas pressure nc^ of the LBM. We verified that the Mach number Ma = \u\max/cs, 
where |u|maa: is the maximum value of the fiuid velocity during evolution, stays 
always smaller than 0.1 so that the fluid results practically incompressible with 
n~ 1. 

The convection-diffusion (12. 5p and the temperature equations are solved by using 
a finite-difference scheme, which is described with details in the paper by Tiribocchi 
et al. (2009). The function (p{r, t) is defined on the same sites of the LBM with the 
same space and time steps. We update (p using an explicit Euler scheme (Strikwerda 
1989). The spatial differential operators are calculated by a second-order finite- 
difference scheme. 

Periodic boundary conditions are used along the a;-direction (and y-direction in 
d = 3) and fiat walls are placed at the lower and upper rows of the lattice along 
the y-direction (z-direction in d = 3), where the temperature is kept at fixed value 
Ttu. Moreover, neutral wetting condition for the order parameter ip is implemented 
at walls imposing that a • ^ip\waii = 0, being a an inward normal unit vector to 
the boundaries. To guarantee exact conservation of the order parameter ip we also 
require that a-\7 {\/^p)\i^aU = Si-V[{V'^)'^p]\waii = 0. The bounce-back rule (Lavallee 
et al. 1991, Desplat et al. 2001) for the LBM is adopted at walls supplemented by 
the prescription for density exact conservation introduced by Lamura & Gonnella 
(2001). 



3. Results and discussion 

In this section we will present the results of our simulations of equations (|2.3p - 
We used the set of parameters b = 0.1, k = -c = -0.03, = 0.005, and 
F = 0.1 corresponding to the lamellar phase in the phase diagram. Lattices had sizes 
256 X 256 in d = 2 and 32 x 32 x 256 in d = 3. The system, initially, is in disordered 
symmetric states with ip fluctuating around zero in the interval [—0.01,0.01]. The 
initial temperature is above with the walls at = 0.9Tc. 

We first show in figure 1 typical configurations of systems with very high vis- 
cosity. In this case we neglected the coupling with the velocity field by solving 
equation (|2.5p without the convective term. Lamellar domains start to form close 
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Figure 1. Configurations of the concentration field ip without hydrodynamics 
with thermal conductivity ^ — 10^^, lO"'^, +00 (from left to right) at times 
t = 55.7 X 10^; 11.7 x 10^; 1.95 x 10^, respectively. 




Figure 2. Configurations of the concentration field (fi at low viscosity (rj = 1.5) 
with thermal conductivity ^ = 10^'', 10^^, +00 (from left to right) at times 
t=UAx 10*; 35.6 x 10*; 2.1 x 10*, respectively. 

to the walls, where the temperature has become lower than the critical value, with 
interfaces perpendicular to the walls in agreement with the imposed neutral wet- 
ting conditions. Then, following the temperature fronts, the regions with phase 
separated domains will extend towards the middle of the system. The pattern mor- 
phology, however, depends on the value of the thermal conductivity or, equivalently, 
on the speed of the temperature fronts. At very large ^, the system remains frozen 
in intertwined noodle-like configurations without long range order, like that on the 
right of figure 1, as observed in other simulations (see, e.g., Gonnella et al. 1997). 
However, by decreasing the value of ^, in the interval [0.1, 1], in spite of the adopted 
boundary conditions, the interfaces tend to follow the temperature fronts, parallel 
to the walls, and exhibit only few defects. The two fronts separating the regions 
with ordered lamellae from the central disordered region will approach each other 
joining at the later times of the simulation. By further decreasing ^, interfaces show 
the tendency to keep the perpendicular orientation with respect to the walls. This 
phenomenology has to be compared with that observed in simple binary mixtures 
where at large £, usual isotropic phase separation occurs (of course, in binary mix- 
tures, there is no frustration and the average size of domains grows by power law, 
see, e.g., Bray 1994), while by lowering the values of the thermal conductivity one 
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Figure 3. Interface configurations in lamellar ordering without hydrodynamics 
with thermal conductivity ^ — 10~^, 1, 10^ (from top to bottom) at times 
t = 78 X 10^; 15 X 10^; 1.1 x 10^, respectively. 



observes domains first parallel and then perpendicular to the cold walls (Tiribocchi 
et al. 2010). 

Going to a lower viscosity {t — b^r] — 1.5, see figure 2), at very large ^, one finds 
that hydrodynamics favours the disentanglement of the network pattern (Gonnella 
et al. 1997) but one still observes extended defects between grains of different ori- 
ented lamellae (see e.g. the grain with vertical orientation on the top-left region of 
the system opposed to the large horizontal underlying domain). One could see in 
systems larger than those here considered that these extended defects would slow 
down the ordering process at very late times (Xu et al. 2005). At lower thermal 
conductivities, in the range [5 x 10"^, 5 x 10~^], one again finds that lamellae are 
well aligned with the thermal fronts. The velocity helps this ordering and, indeed, 
very ordered defect-free lamellae can be observed (compare the central snapshots 
of figure 1 and figure 2) in the regions close to the walls. In the central part, even 
if structures can be observed, the local values of the field <p are nor in equilibrium, 
and further evolution will show a completely ordered lamellar state. Finally, by 
further decreasing the value of ^, one sees in the left snapshot of figure 2 that hy- 
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dro dynamics clearly helps the tendency of the system to prefer the perpendicular 
orientation. 

Wc also studied three-dimensional lamellar systems without coupling to the 
velocity field. At large values of ^ the system evolves towards metastable configu- 
rations ordered only on short scales, similarly to what occurs in two-dimensional 
systems, as it can be seen in the bottom snapshot of figure 3. At lower values of 
^ (0.1 < ^ < 10) one can see stacks of lamellae parallel to the thermal fronts with 
very few defects (central snapshot of figure 3), while, again similarly to the two- 
dimensional case of figure 1, perpendicular order is lost far from the walls at very 
low values of ^. 

4. Conclusions 

In this work we have studied the ordering of a lamellar phase where the temperature 
cold fronts diffusively move from external walls. Phase separation therefore starts 
close to the walls but develops in the middle of the systems with a phenomenology 
depending on the speed of the thermal fronts. In an intermediate range of values of 
thermal conductivity ^, lamellae appear ordered and parallel to the thermal fronts 
while a perpendicular orientation is preferred at very low values of ^. Hydrody- 
namics favours these preferred ordering and one can conclude that quenching with 
moving temperature fronts are effective in producing lamellar states with very few 
defects. A full analysis of the three-dimensional case will complete this study. 

We acknowledge support from CASPUR for "HPC Grant 2010". 
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